Generation of dark-bright soliton trains in superfluid-superfluid counterflow 
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The dynamics of two penetrating superfluids exhibit an intriguing variety of nonhnear effects. 
Using two distinguishable components of a Bose-Einstein condensate, we investigate the counterflow 
of two superfluids in a narrow channel. We present the first experimental observation of trains of 
dark-bright solitons generated by the counterflow. Our observations are theoretically interpreted 
by three-dimensional numerical simulations for the coupled Gross-Pitaevskii (GP) equations and 
the analysis of a jump in the two relatively flowing components' densities. Counterflow induced 
modulational instability for this miscible system is identified as the central process in the dynamics. 
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Nonlinear structures in dilute-gas Bose-Einstein con- 
densates (BECs) have been the focus of intense research 
efforts, deepening our understanding of quantum dy- 
namics and providing intriguing parallels between atomic 
physics, condensed matter and optical systems. For su- 
perfluids that are confined in a narrow channel, one of 
the most prominent phenomena of nonlinear behavior is 
the existence of solitons in which a tendency to disperse 
is counterbalanced by the nonlinearities of the system. In 
single-component BECs, dark and bright solitons, form- 
ing local density suppressions and local bumps in the 
density, resp., have attracted great interest [l|. In two- 
component BECs, the dynamics are even richer as a new 
degree of freedom, the relative flow between the two com- 
ponents, is possible. 

In this Letter, we investigate novel dynamics of 
superfluid-superfluid counterflow, which is in contrast 
to the extensively studied counterflow of a superfluid 
and normal fluid in liquid helium [2|. Previous theo- 
retical analysis has demonstrated that spatially uniform, 
counterfiowing superfluids exhibit modulational instabil- 
ity (MI) when the relative speed exceeds a critical value 
[3|. Modulational instability is characterized by a rapid 
growth of long wavelength, small amplitude perturba- 
tions to a carrier wave into large amplitude modulations. 
The growth is due to the nonlinearity in the system \^. 
Our experiments and analysis reveal that by carefully 
tuning the relative speed slightly above the critical value, 
we can enhance large amplitude density modulations at 
the overlap interface between two nonlinearly coupled 
BEC components while mitigating the effects of MI in 
the slowly varying background regions. A dark-bright 
soliton train then results. 

In two previous experiments, individual dark-bright 
solitons were engineered in two stationary components 
using a wavefunction engineering technique J5|, |6|. In 
our experiment we find that trains of dark-bright solitons 
can occur quite naturally in superfluid counterflow. This 
novel method of generating dark-bright solitons turns 



out to be robust and repeatable. In single-component, 
attractive BECs the formation of a bright soliton train 
from an initial density jump has been predicted [7|. How- 
ever, both condensate collapse and the effects of MI in 
the density background must be avoided, placing restric- 
tions on the confinement geometry and diluteness of the 
single-component condensate. In contrast, the proper- 
ties of counterflow in miscible, two-component BECs, as 
we show, enable the observation of trains consisting of 
ten or more dark-bright solitons in BECs with a large 
number of atoms. We also note that modulated soliton 
trains have been studied extensively in single-component, 
modulationally stable repulsive BECs where supersonic 
flow supports the generation of dispersive shock waves 
[8J. The dark-bright soliton train we study in this work 
occurs when one of the system's sound speeds becomes 
complex so that the standard definition of supersonic flow 
does not apply. 

Our experiments are conducted with BECs confined 
in a single-beam optical dipole trap [9]. We start with 
an initially perfectly overlapped mixture of atoms in the 
|F, m^) = |1,1) and |2,2) hyperfine states of ^^Rb, with 
a total of about 450000 atoms. The scattering lengths for 
the two states used in our experiment are estimated to be 
ail = 100.40 a.u. and a22 ^ <^i2 = 98.98 a.u. [lO|. Here 
ail and a22 denote the single species scattering length 
for the |1, 1) and |2,2) state, respectively, and ai2 is the 
interspecies scattering length. Mean field theory predicts 
that a mixture is miscible if ai2 < y'an • a22 Ill4l3|. 
Therefore our system is predicted to be weakly miscible. 
In contrast, previous studies of two-component binary 
^^Rb BECs concentrated mostly on the states which are 
immiscible IJ, ll5| , with the notable exception of Weld 
et al. [Tel. 

When the overlapped mixture is allowed to evolve in 
the trap, we observe no phase separation over the experi- 
mental timescale of several seconds. This is in agreement 
with the predicted miscibility of the two components and 
is demonstrated in Fig. [TJa-c). The upper cloud of each 
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FIG. 1: Time evolution of an initial perfectly overlapped mix- 
ture without (a-c) and with (d) an applied axial magnetic 
gradient. Images taken after (a) 100 ms, (b) 1 sec and (c, 
d) 9 sec of in-trap evolution. 



image throughout this work shows the atoms in the |2, 2) 
state at a time 7 ms after a sudden turn-off of the opti- 
cal trap (and, where applicable, of any applied magnetic 
gradients), while the lower cloud, taken during the same 
experimental run, shows the atoms in the |1, 1) state af- 
ter 8 ms of expansion [9y . During their in-trap evolution, 
these clouds are overlapped in the vertical direction. The 
dominant effect of the time evolution in Fig. [T]^a-c) is a 
slow decay of the atom number over time. For single 
component BECs, we have measured an exponential BEC 
lifetime of over 50 sec for the |1, 1) state and 14 sec for the 
1 2, 2) state in our dipole trap. Motion induced by changes 
of mean field pressure during the decay may be responsi- 
ble for a small scale roughness of both components which 
becomes visible after several seconds (Fig. [Th). 

The situation changes when a small magnetic gradient 
is applied along the long axis of the trap. Due to Zeeman 
shifts, the gradient leads to a force in opposite directions 
for each component, or equivalently to a differential shift 
between the harmonic potentials along the long axis of 
the trap. This causes the two components to accelerate 
in opposite directions and induces counterflow. In all 
images where a magnetic gradient is applied, the gradient 
is chosen such that the |2,2) state is pulled to the right 
and the |1, 1) to the left. An example is shown in Fig.[l]i 
where a gradient leading to a calculated differential trap 
shift of 60 fim. was applied for 9 sec, leading to nearly 
complete demixing of the two components. 

In the following we investigate the dynamics induced 
by small gradients and show how they can be exploited 
to create dark-bright soliton trains. In Fig. [2] an ini- 
tially overlapped mixture of 30% of the atoms in the |2, 2) 
state and 70% in the |1,1) state is used. A small mag- 
netic gradient in the axial direction is linearly ramped 
on over a timescale of 1 sec, leading to a calculated trap 
separation for the two species of only about three mi- 
crons. After the end of this ramp, the gradient is held 
constant. In the subsequent evolution, individual stripes 
break off from the left edge of the |2, 2) component, and 
perfectly aligned dark notches appear in the |1, 1) com- 
ponent (Fig. EJa)). The predominantly uniform widths 
of the observed stripes and notches, their long lifetime 
of several seconds in the absence of a magnetic gradient. 



as well as their dynamics resembling individual stable 
entities (see Fig. |4] and below) are strong experimental 
indications that the observed features are indeed dark- 
bright solitons. By reducing the initial number of atoms 
in the component forming the bright soliton, we have also 
been able to rehably produce one individual dark-bright 
soliton and observe its oscillation in trap [17|, similar to 
the dynamics observed in [6|. 

The observed soliton formation is reproduced by three- 
dimensional (3D) numerical simulations of the two- 
component Gross-Pitaevskii (GP) equations (Fig. [2]^b- 
e)) [9|. Parameters used for the GP equations are the 
experimental values. These values lead to dynamics that 
closely match the experiment, as shown in Fig. [2fa-c) 
with a moderate time delay. Our numerical calculations 
suggest that the time delay may be due to uncertainties in 
the estimated magnetic field gradient induced trap shifts. 
The experimentally invoked free expansion directly be- 
fore imaging the condensate was not performed in the 
numerical simulations. 

Numerical results for the quantum mechanical phases 
of the two wavefunctions describing the components are 
shown in Fig. [2je). The nearly linear phase behavior on 
the right (at x > 50jj.m) indicates a smooth counterflow 
of the two components. In the soliton region, the phase 
jumps across the dark solitons as well as the phase gra- 
dients in the bright component vary slightly, so that the 
dark-bright solitons are moving relative to one another 
which eventually leads to dark-bright soliton interactions, 
see [9]. 

The soliton train formation can be qualitatively un- 
derstood by appealing to the hydrodynamic formulation 
of the mean- field, coupled GP equations in (1+1) dimen- 
sions 
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here given in non-dimensional form with aj = ai2/ajj, pj 
and Uj^ j = 1,2 the density and phase gradient (super- 
fluid velocity) of the j^^ component, respectively. Equa- 
tion ([T|) models the dynamics of a highly elongated cigar 
shaped trap {ujx ^ ojy ^ ujz where ujx ioOy) is the trans- 
verse trap frequency in the horizontal (vertical) plane 
and ujz is the axial trap frequency) with axial confine- 
ment neglected [l| . Distance is in units of the transverse 
harmonic oscillator length y/h/{mujx) (rn is the particle 
mass). Time is in units of l/uox and the 3D densities are 
approximated by the harmonic oscillator ground state via 
p^-(2;,t)exp(-x2 - ^y)/{2iiajjal). 

By considering small perturbations proportional to 
^i{Kz-ujt) £q^ uniform counterflow with densities pj and 
velocities ui = —v/2, U2 = v/2, Ref. [3[ demonstrated 
modulational instability {Im. uj{k) > 0) for v larger than a 
critical velocity 'i;cr with a maximum growth rate Im cjmax 



(a) 
1000 ms 

1100 ms 

1200 ms 

1300 ms 

1400 ms 

1500 msl 



(b) 








1052 ms 



1300 ms 




c 



(C) 



ZIK^ 



-300 




x[^m] 



H-300 




100 



100 



x[nnn] 



FIG. 2: Formation of dark-bright soliton trains during superfluid-superfluid counterflow. A mixture with 70% of all atoms in 
the |1,1) state (lower cloud in (a) and (c); black line in (b),(d),(e)) and 30% in the |2,2) state (upper cloud and red lines, 
resp.) is created. An axial gradient is ramped on over 1 sec and held constant thereafter, (a) Experimental images of soliton 
train formation. Times are measured from start of ramp, (b) Integrated cross sections of 3D numerical simulations, revealing a 
gradually steepening overlap interface and subsequent soliton formation, (c) Integrated density plot of 3D numerics at 1751 ms. 
Field of view 627 /xm x 17 /im. (d) Zoomed-in view of soliton train in (b) at 1751 ms. The blue line shows the total density. 
(e) Phase behavior for region shown in (d). 



and associated wavenumber A^max- We have repeated the 
calculation and find the additional result 



\/pl{l - (JlCr2) < ^cr < 2y Pl(l - y/(JlCF2), pi > p2 , 

(2) 
the lower bound being valid for small P2/P1 and the upper 
bound applicable for p2 ^ pi. The scattering lengths of 
the binary system considered here give 0.1 19y^ < 'Ucr < 
0.168y^. Typical densities for the experiments in Fig. [2] 
give pi -\- p2 = 4.3, P2/P1 = 0.3 leading to Va = 0.25 
{^ 0.22 mm/s). 

Figure [2] shows that the dark-bright soliton train forms 
at the overlap interface of the two components while 
approximately maintaining constant total density. We 
model this by numerically solving eq. ([T]) for an initial 
jump in density that maintains pi -\- p2 = 4.3 (dotted 
curves in[3l^a,b)) with a uniform counterflow: Ui = —v/2 
and U2 = v/2. For subcritical cases < v < ^cr, the 
evolution consists of an expanding rarefaction wave with 
weak oscillations on the right edge (Fig. [3l^a); solid line 
for '^ = 0, dashed line for v = 0.17) and correspond- 
ing scaled relative speeds \ui — U2\l ^fpi (Fig. [3l^c) solid, 
dashed) below critical (in Fig. [3l^c,d), the bounds (|2]) on 
the critical velocities are indicated by the dotted lines). 
When the initial relative speed is supercritical, a dark- 
bright soliton train forms at the initial jump (Fig. [3fb), 
V = 0.32). The relative speed within some regions of 
the soliton train significantly exceeds ^cr as shown in 
Fig. [3l^d) suggesting that counterflow induced MI has the 
effect of enhancing soliton formation. Because the initial 



relative speed v was taken just slightly above v^^^ the 
maximum growth rate ImcJmax = 0.0077 and associated 
wavenumber /^max = 0.13 for unstable perturbations to 
the uniform state in the far field are small (Figs. [3l^e,f)). 
Therefore, MI in the background counterflow far from 
the jump does not develop appreciable magnitude over 
the timescale of soliton train formation, in contrast to 
the dynamics with v ^ ^cr that we investigate in [18|. 

This MI assisted soliton formation technique allows us 
to create dark-bright solitons in a well-controlled and re- 
peat able manner, as is evidenced by the fact that all 
images of Fig. [2^ form a very consistent sequence even 
though they were taken during different runs of the ex- 
periment. In addition to repeatability, future studies may 
also require a long lifetime of the solitons. In single com- 
ponent BECs, achieving long lifetimes of dark solitons 
has proven difficult as they are subject to a transverse 
instability [5|, |8|. Only recently have dark soliton life- 
times of up to 2.8 sec been achieved [6[. It has been con- 
jectured [l9[ and numerically conffimed [20| that dark- 
bright solitons are more stable to transverse perturba- 
tions than dark solitons. Experimentally, we indeed ob- 
serve long lifetimes of several seconds for the dark-bright 
solitons after the magnetic gradient is turned off. The 
solitons act as individual entities and can move through 
the BEG, maintaining their shape for a relatively long 
time. We demonstrate this by starting from a situation 
as in Fig. [2fa) at 1.5 sec, where a train of solitons has 
been created after the application of an axial magnetic 
gradient. When the gradient is subsequently turned off. 
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FIG. 3: densities (a,b) and relative speeds (c,d) at t = 462. 
(e,f): maximum growth rate and associated wavenumber of 
unstable perturbations. Vertical lines correspond to subcrit- 
ical (dotted) and supercritical (dashed) cases of (a-d). See 
text for further details. 
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FIG. 4: Motion of dark-bright solitons. After creating dark- 
bright solitons as in Fig. [2l^a), the axial magnetic gradient is 
removed and the mixture is allowed to evolve in trap for the 
evolution times given next to each image. The solitons spread 
out through the mixture. 
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the dark-bright solitons move through the BEG while 
approximately maintaining their narrow widths (Fig. |4]). 
The bright and dark part of each individual soliton re- 
main aligned relative to each other, but any regularity in 
the spacing between solitons is lost. The number of visi- 
ble solitons decreases over time, but even after 2.5 sec sev- 
eral solitons are still visible, as in Fig. HJ^d). Simulations 
[9[ suggest that soliton interactions may be the cause of 
this decay. In Fig. HJ^c), little diffuse cloudlets of atoms 
in the |2,2) state are visible in addition to some soli- 
tons, and corresponding small suppressions of the density 
in the |1,1) components can be detected. We interpret 
these features as the decay products of dark-bright soli- 
tons, marking the end of their life cycle. 

In conclusion, we have observed dark-bright soliton 
trains in the counterfiow of two miscible superfluids. The 
soliton train is formed due to relative motion above the 
critical value for modulational instability. By inducing 
relative speeds slightly above critical, we can avoid the 
onset of MI throughout the superfluids over the time 
scales of soliton train formation. Together with the long 
lifetime of the observed dark-bright solitons, this opens 
the door to future experiments with these interesting co- 
herent nonlinear structures. While the dynamics con- 
sidered in this Letter are effectively one-dimensional, a 
very recent theoretical analysis has shown that super- 
fluid counterfiow in higher dimensions can lead to bi- 
nary quantum turbulence, providing another example of 
the exceptional dynamical richness of the two-component 
system j21l] . 
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